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A ID, steady-state model for Brillouin and Raman backscatter from an inhomogeneous plasma 
is presented. The daughter plasma waves are treated in the strong damping limit, and have ampli- 
tudes given by the (linear) kinetic response to the ponderomotive drive. Pump depletion, inverse- 
bremsstrahlung damping, bremsstrahlung emission, Thomson scattering off density fluctuations, and 
whole-beam focusing are included. The numerical code deplete, which implements this model, is 
described. The model is compared with traditional linear gain calculations, as well as "plane-wave" 
simulations with the paraxial propagation code pF3D. Comparisons with Brillouin-scattering exper- 
iments at the OMEGA Laser Facility [T. R. Boehly et al., Opt. Commun. 133, p. 495 (1997)] show 
that laser speckles greatly enhance the reflectivity over the deplete results. An approximate upper 
bound on this enhancement, motivated by phase conjugation, is given by doubling the deplete 
coupling coefficient. Analysis with deplete of an ignition design for the National Ignition Facility 
(NIF) [J. A. Paisner, E. M. Campbell, and W. J. Hogan, Fusion Technol. 26, p. 755 (1994)], with 
a peak radiation temperature of 285 eV, shows encouragingly low reflectivity. Re-absorption of 
Raman light is seen to be significant in this design. 
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I. INTRODUCTION 

Laser-plasma interaction (LPI) [T] is an important 
plasma-physics problem which poses serious challenges to 
theoretical modeling. LPI is the basis of several applica- 
tions, including laser-based particle acceleration [3] and 
the backward Raman amplifier [5] . Moreover, for inertial 
confinement fusion (ICF)[H|5] to succeed, LPI must not 
be so active that it prevents the desired laser energy from 
being delivered to the target, with the desired spatial and 
temporal behavior. This paper focuses on modeling the 
backscatter instabilities, where a laser light wave (mode 
0) decays into a backscattered light wave (mode 1) and 
a plasma wave (mode 2) . In stimulated Raman scatter- 
ing (SRS) and stimulated Brillouin scattering (SBS), the 
plasma wave is, respectively, an electron plasma wave 
and an ion acoustic wave. These LPI processes pose a 
serious risk to indirect-drive ICF [5]. 

A wide array of computational tools is used to model 
LPI, ranging from rapid (~secs) calculations of lin- 
ear gains along ID profiles to massively-parallel kinetic 
particle-in-cell simulations. We present here a new tool, 
called deplete, to the less computationally expensive 
end of this spectrum, deplete solves for the pump inten- 
sity and scattered-wave spectral density for a set of scat- 
tered frequencies, in steady-state, along a ID profile of 
plasma conditions. Pump depletion is included, and the 
plasma waves are assumed to be in the strong damping 
limit (i.e., they do not advect). Fully kinetic (although 
linear) formulas are used for various quantities like the 
coupling coefficient. Bremsstrahlung noise and damping, 



as well as Thomson scattering (TS), are included. The 
deplete model, especially the noise sources, in some 
ways resembles that of Ref. [B]. Other similar works 
which have influenced our thinking, and use ID coupled- 
mode equations, are Refs. [7]- [8]. 

deplete is a ID model, but the plasma conditions are 
generally found by tracing 3D geometric-optics ray paths 
through the output of a radiation-hydrodynamics code. 
We therefore call this combined approach to studying 
LPI a ray-based one. Details of this methodology, and 
its limits, are discussed in Sec. |III| 

deplete is similar to the code newlip, which cal- 
culates linear gains for SRS and SBS along ID profiles 
(newlip is discussed here in Appendix A). Both codes 
take seconds to analyze one profile from the laser en- 
trance to the high-Z wall in an ICF ignition design. 
However, deplete includes substantially more physics 
than newlip, such as pump depletion, noise sources, 
and re-absorption of scattered light, deplete more- 
over provides pump and scattered intensities, which un- 
like gains can be directly compared with experiment and 
more sophisticated LPI codes. Despite its simplicity, de- 
plete agrees well in certain cases with results from the 
3D paraxial laser propagation code pF3D. This is quite 
promising given DEPLETE's much lower computing cost. 

There is important physics which deplete does not 
capture, with laser speckles or hot spots being one of 
the most important. Recent SBS experiments [21 EH] at 
the OMEGA Laser Facility [TT] show good agreement be- 
tween measured reflectivity and pF3D predictions, while 
deplete gives a lower value. This is due to the speckle 
pattern of the phase plate smoothed lasers. Sec. VIII 
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describes one approximate way to bound the speckle en- 
hancement by doubling the coupling coefficient; the re- 
sulting deplete reflectivity always exceeds the exper- 
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imental level. A more sophisticated idea for handling 
speckles is outlined in the conclusion. Additional beam 
smoothing, like polarization smoothing (PS) and smooth- 
ing by spectral dispersion (SSD), reduce the effective 
speckle intensity and can reduce the reflectivity even be- 
low the speckle-free deplete level. 

The paper is organized as follows. Section |TT] de- 
rives the governing equations for the pump intensity and 
scattered-wave spectral density. Our ray-based method- 
ology and model limits are discussed in Sec. The 
numerical method is given in Sec. |IV[ including a quasi- 
analytic solution for the coupling- Thomson step. Section 
[V]compares DEPLETE with NEWLIP linear gains and pF3D 
"plane-wave" simulations on prescribed profiles. The re- 
lationship between Thomson scattering and linear gain 
is discussed in Sec 



where i = 0, 1 for the pump and scattered wave, respec- 
tively. Ai is the slowly- varying complex envelope, and 
we use the dimensionless a, = eAi/m e c. tpi(z,t) is the 
rapidly-varying phase with k% = d z ipi and u>i = —dtipi- 
Let <7 2 = ki/\ki\ with <jq = 02 = +1 and a\ = — 1 (ap- 
propriate for backscatter) . Thermal fluctuations give rise 
to both light waves and plasma waves. However, upon 
appropriate averaging the field amplitudes of these fluc- 
tuations vanish (but their mean squares do not). The 
amplitudes Ai (and nj% below) represent only the coher- 
ent, and not the noise, components of the fields. We 
insert a bremsstrahlung noise source and TS to the in- 
tensity equations below. 

From the Maxwell equations, and conservation of 
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canonical transverse momentum m e v ye 



eA y , we find 



to the experimental and pF3D SBS reflectivities in recent 
OMEGA shots. Sec. |VIII| presents deplete analysis of 
an ignition design with a 285 eV radiation temperature 
for the National Ignition Facility (NIF) [12]. In particu- 
lar, we show the effect of scattered light re-absorption and 
put a bound on speckle enhancement. We conclude and 
discuss future prospects in Sec. |IX| A review of NEWLIP 
and its linear gain is presented in Appendix A. Appendix 
B details the numerics of DEPLETE's coupling- Thomson 
step. 



II. GOVERNING EQUATIONS 

We derive coupled- mode equations, in time and one 
space dimension, for the slowly-varying wave envelopes, 
and find the resulting intensity equations. We do this for 
the light waves first, and then the plasma wave in the 
strong damping limit. Since our approach is standard we 
summarize some steps. We take these equations in steady 
state to apply independently at each scattered frequency, 
and transition to a spectrum of scattered light per an- 
gular frequency. This may be viewed as a "completely 
incoherent" treatment of the scattered light at different 
frequencies. Bremsstrahlung damping and fluctuations, 
and TS, are then added phenomenologically. Focusing of 
the whole beam is finally accounted for, giving the system 
deplete solves. This section culminates in the deplete 
system, Eqs. ( 54|55 l, on which some readers may wish to 
focus. 



we compare deplete A y = (A a + A t ) ■ y satisfies 



A. light-wave action equations 

Let z be distance along the profile, and assume all 
wave vectors and gradients are in z (d x — d y = 0). 
z = is taken as the left edge of the domain (the 
"laser entrance" ) , where we specify the right-moving 
pump laser; we also specify boundary values for the left- 
moving backscattered wave at the right edge z = L z . The 
light waves are linearly polarized in y and represented 
by their vector potentials A; = (l/2)A;(z, t)ye 1 ^ i + cc, 
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(1) 



hj = rij + Nj2 is the total number density for species 
j (j = e for electrons, i for an ion species), Nj2 — 
(l/2)nj2e 1 ^ 2 + cc, and nj2 is the slowly- varying plasma- 
wave envelope. We define u> p j = [njZ 2 e 2 / 'e^raA 1 / 2 , 

VTj = [Tj/rrij] 1 / 2 and Xnj = ^Tj/^pj, with Zj the charge 
state. As usual, the massive ions are treated as fixed in 
the transverse current. (We look forward to a circum- 
stance where a positively-charged species must be con- 
sidered mobile, such as an electron-positron plasma!) 

Following, e.g., Ref. [13], we introduce the small pa- 
rameter S ~ ui~ 1 dtlnX ~ k~ 1 d x lnX for X = Ai,ki, 
etc. We order d t ,d x ~ S, ipi ~ S^ 1 , and the right-hand 
side of Eq. ([TJ ~ d. To order 5°, we obtain the free- wave 
dispersion relation 
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For the steady-state conditions considered below we take 
u>i to be constant and find the eikonal ck^x) = aiijiUJi 
with rji = \V~n e jn C i[ L l 2 and n C i = uj 2 eom e /e 2 the critical 
density of mode i. Also, the group velocity is v g i = atrjiC. 

Assuming perfect phase matching (fco = k% + k-2, luq = 
uii + cj 2 ), the resonant order S terms in Eq. ([!]) yield the 
envelope equations: 



^0^0 — — 



4 w n e 
i Up e n%2 
4 u>i n e 



a . 



(3) 
(4) 



The operator Li = d t + v g id z + (l/2w i )(9 t a;i + c 2 d z ki). 
Our quasi-monochromatic light waves (i = 0, 1) have 
action density [T4j iVj = (m e /87rr e )u;iOia* where r e = 
e 2 /4:TT€oin e c 2 « 2.82 fm. We also define the (positive) 
action flux Zj = 7Vj|?; g j| and intensity = WjZj. In prac- 
tical units, 



h\ 2 



(5) 
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where A, = 2Trc/u>i and P em = (w/2)m e c 3 /r e « 1.37 x 
10 18 W-cm -2 • ^m 2 . We form Eq. |3])xa|$ + cc and 
Eq. ^xal+cc to find 



d t N 



d z Z 
J 



1 



d z Z x = J 



m e c Im[aj5airo e 2]. 



(6) 
(7) 



B. plasma- wave action equations 



We describe the plasma waves following the dielectric 
operator approach of Cohen and Kaufman |15j : 



e(ui 2 + id t , k 2 - id z )n 2 = 



^pnd i 



(8) 



« P nd = Xe{^k2)y^-n e a al.{^) 

pe 

The charge-density fluctuation n 2 = ~ n e2 + ^2i ZiUi 2 ex- 
periences a ponderomotive drive n pn d- w 2 = lo 2 — k 2 ■ u is 
the Doppler-shifted plasma-wave frequency in the frame 
of the plasma flow u (u> 2 is in the lab frame) . e = 1 + x is 
an operator, where the time and space derivatives reflect 
envelope evolution and x = J2j Xj 1S the total suscepti- 
bility. Xe in n pnd is simply a function, not an operator. 
Xj is the (linear) kinetic, collisionless susceptibility of 
Maxwellian species j: 
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e * erfc(— iQ is the plasma dispersion func- 



tion |16j and erfc is the complimentary error function |17j . 
Gauss's law relates n 2 and rij 2 : 



n e2 = 
n l2 = 



-(1 + Xi)n>2, 

1 m e e 

Zi m,i Xe 



-Xi 



"2 



Xt 



(11) 

(12) 
(13) 



with xi = T^iXi- For SRS, where the ion motion is 
negligible, we usually take 1 + xi — * 1 to save computing 
time. 

Expanding e for slow envelope variation, and retaining 
only e r = Ree in the derivatives, gives 



[dt + v g2 d z + + iSuj 2 ] n 2 



^pnd 



(14) 



e = de r /duj 2 , e' = de r /dk 2 , v g2 = —e'/e is the plasma- 
wave group velocity, 1S2 = Im[e]/e is the damping rate, 
and Suj2 = — e r /e is the phase detuning. 

We now assume the plasma wave is in the strong damp- 
ing limit, where its advection is neglected: \v g 2d z ri2\ <C 
\v 2 + i<5cij2 1 1^2 1 - This implies the instability is below its 
absolute threshold so that steady-state solutions are ac- 
cessible. Also going to steady-state, we find 



e(uj 2 ,k 2 )n 2 = n pnd . 



(15) 



Replacing n e 2 via Eqs. (11 1 and (15 1 yields 



J — ojqTiZoZi 
The coupling coefficient r\ is 



9t 



T s lm 



<l + Xl) 



Ts.9r 



2nr e 1 k 2 



m e c 2 u k \ki\ 

|1 + X/| 2lm Xe + IXe| 2 ImX/. 



(16) 



(17) 

(18) 
(19) 



The second form of f i exhibits the resonance for |e| <C 1. 
The over-tilde on Ti indicates it will be modified below 
to account for beam focusing. T±, and thus J, are usually 
positive. We now have a closed system for modes and 
1, with no independent equation for mode 2: 

d t N + d z Z = -uqTiZqZx, (20) 
dtNt-dzZ! = wofi-ZbZi. (21) 



C. Steady-state equations for a spectrum of 
scattered waves 

We transition to steady state (dt = 0) and work with 
intensities. Since we have assumed d z uii = 0, we multiply 



Eq. (20) by ojq and Eq. (21 1 by lo\ to obtain 



d z Io = 
-d z h = Tiloh 



W °f T T 



(22) 
(23) 



Here and elsewhere, d x f(x) denotes the ordinary deriva- 
tive of a function of one variable, while d x f denotes the 
partial derivative of a function of several variables. 

The bremsstrahlung source and TS are expressed in 
terms of spectral density ii(z,wi) (intensity per angular 
frequency). The scattered intensity is then I\ — J du>\ i\. 
We take Eq. ( 23 ) to apply independently at each u>i , and 



integrate the coupling term in Eq. (22 1, to find 



d z I 



duii — Til ii, 



-d z i x = Txloh. 



(24) 
(25) 



This is a totally incoherent treatment of the scattered 
light at different frequencies, and is unrealistic to the 
extent there is spectral "leakage" between nearby lo\ in- 
tervals due to, e.g., envelope evolution. 



D. Bremsstrahlung source and damping 

We incorporate electron-ion inverse-bremsstrahlung 
light- wave damping (kq and ki) phenomenologically for 
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modes and 1, as well as bremsstrahlung noise (Si) for 
mode 1, to find 



d z I = -kqIq- I duJi— fi/oil) 

UJi 



-d z ii 



-K\ii + Ei + Til ii. 



(26) 
(27) 



As for f i, the over-tilde on Si denotes it will be modified 
due to focusing. 

Iq and i\ represent integrals over solid angles in k 
space, which we now specify. Absolute solid angles are 
needed in the noise sources, and cannot be simply scaled 
away, because scattered intensities determine pump de- 
pletion. We follow closely Bekefi's book [18] in this sec- 
tion. We take Jj = fijl^n for i = 0, 1 (see Sees. 1.6 and 
1.7 of Bekefi). I^q is the intensity per solid angle interval 
dQ, in k space, which we assume is constant over the solid 
angle f2j that participates in the scattering. f2j is the lo- 
cal (in z) solid angle in the plasma, which we express in 
terms of a cone half-angle 9 p ^ as 

fli = 2%{l- cos 6 Pti ). (28) 

From Snell's law, 9 P ^ varies with z according to 



cos 9 V 




cos ( 
r 2 sin 2 ^lV2 



otherwise. 



(29) 



n c i cos 2 9 V is the "critical density" above which we cut off 
backscatter (l\ = Si = m = 0). 9 V is a "vacuum" cone 
angle, which we find from the solid angle in the beam's F- 
cone (for simplicity we use the same solid angle for pump 
and scattered light). This is reasonable if the scattering 
mostly occurs in laser speckles that are near diffraction- 
limited. In terms of laser optics F-number F, 



cos 9 V 



i 



i 

4F2 



-1/2 



2n(l — cosd v ) 



s 1 - 

7T 



1 

8F 2 "' 



(30) 
(31) 



The approximate forms apply for F 3> 1. 

The upshot of the solid angle discussion (see especially 
Eq. (1.133) of Bekefi) is 



Si = Qij(u)i), 



(32) 



where j(ui) is the emission coefficient, per dfl and in one 
polarization (see p. 134 of Bekefi): 



u} pe m e r e 



(33) 



127r 3 V27TWTe < 

j £ ions 

In A e j is sometimes called the Gaunt factor and resembles 
the Coulomb logarithm, although it arises in calculations 
without ad hoc cutoffs on impact parameter integrals (see 
Chap. 3 of Bekefi). For the case a;, > uj pe , Bekefi finds 
) where 



A e j = VTe/^LO.b 



h . _ J 4 S^T t 
K -~ \ ft)"'' 



7 . ,,. r,i. t (- 
^ j 1 e j< 



if T e > 77 Z] eV, 
otherwise. 



The first, high-T e case typically applies for hohlraum con- 
ditions. The numerical pre- factors come from a detailed 
binary-collision calculation, and 7 = e c ss 1.781 where 
C « 0.577 is the Euler-Mascheroni constant. Our ex- 
pression for j does not include the enhanced emission for 
uji fa uj pe due to collective effects [19] . 

We find the absorption coefficient via Kirchoff 's law 
(see Bekefi Sec. 2.3): 



(35) 



Our Ki equals Bekefi's a u . B v is the vacuum blackbody 
spectrum for one polarization, with units dI/(du)dVL): 



B v {u) 



87r 3 c 2 e huj/T e 



0J 2 T e . 



(36) 
(37) 



j given above was found for collision durations short com- 
pared to the light-wave period, which entails the Jeans 
limit fioj -C T e . We therefore use the approximate form 
of B v to obtain 



V2 a t 



r e cr]i 



(38) 



jGions 



For an optically thick plasma (d z i\ = 0) with no pump 
(la = 0), we obtain for i\ from Eq. (|27l) the fluctuation 
level 



,-OT. 



•OT 



Si 



f 



(39) 



/ and Si are defined in Sec. II F We thus recover the 
blackbody spectrum, required by Kirchoff's law. The 
factor r/l that usually appears in the blackbody spectrum 
in a plasma is absent due to our treatment of solid angles. 



E. Thomson scattering 

Thomson scattering (TS) refers to scattering off 
plasma- wave fluctuations resulting from particle discrete- 
ness (|20J, p. 308). Had we retained a separate plasma 
wave equation, the fluctuations would appear in it as 
Cerenkov emission [6]. It is an important noise source 
for backscatter, especially for SBS. We express Api, the 
TS scattered power increment per du>i per d£li (k\ solid 
angle), within a thin slab of width Az, as 



da 
duiidCli 



da 



du>idfli 



n e A(z)Azipr' l 



S 

2tt' 



(40) 
(41) 



A(z) is the beam area, defined in Sec. II F ijj = 1 — 
sin 2 9 S sin 2 9 a is a geometric factor. 9 S is the angle be- 
tween kg and R, the vector from source to "observation 
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point" . For a beam with large F, 6 a ~ 9 V -C 1. a is the 
angle between R and the pump polarization. We usually 
take if) = 1. 

The form factor S 1 (units of time) is from Eq. (138) of 
Ref. |20j . valid for arbitrary (non-Maxwellian) distribu- 
tions, generalized to multiple ion species: 



2?r 



S(k,u) = \l + X i\ 2 Fe 



| Xe | 2 Y: ^F 3 (42) 



Fj = / d 3 v fji^Siuj + k -v). 



(43) 



/j is the distribution function of species j (J d 3 v fj = 1). 
For a Maxwellian, 



F< = 



1 



kvTj v2~7r 



e -c? = (^i mx ., 



(44) 



r 5 = ,g T = |1 + x/| 2 lmxe 



^F lul ^ 



and 

^i 2 

2.:/-v\/,, i 

This form agrees with the multiple-ion result in Eq. (3) 
of Ref. [H] • Henceforth we assume Maxwellian distribu- 
tions. 



From Eq. (|40| we form a differential equation for i\ 

f2i Api 



that describes TS: 

<9 2 ii| TS = txIq 



Tl 



AI Az 



(46) 



Since TS transfers energy from the pump to the scattered 
waves, we include it in both equations: 



d z I = -k I - / dwi— I (rx +f 1 i 1 ), 
J wi 

-d z ii = -K 1 i 1 + t 1 +I (T 1 + f 1 i 1 ). 
For conevience we write t\ as 



2 L S(k2,u' 2 ) _ r s g T 



2tt 



t S = n e r e . 



(47) 
(48) 

(49) 
(50) 



Ti is always positive, while rg and g T have the same sign 
as lu' 2 (which can be negative for IAW's when the plasma 
flow is supersonic along ko). 

It is useful to note that i T = Ti/T\ sometimes plays 
the role of an effective seed level for i\. 

■ =IL — Ts9t 
Ti T s gr' 

For the special case — T ei we have g T = <?r and i T is 
independent of \j '■ 



(51) 



TS 



T T e k \ki\, T t = T e 



(52) 



T S (2tt)3 J 2 

This fact is used in Sec. IVII to discuss the relation of TS 
to linear gain. 



F. Whole-beam focusing 

We wish to incorporate the effects of whole-beam fo- 
cusing in a simple way. The equations as written hold 
locally in z, but do not model focusing. To do this, 
we treat the transverse intensity patterns of Iq and 1\ 
to be uniform flattops of varying area A(z). The beam 
focuses at the focal spot zp, where A attains its min- 
imum A(zp)- Let Ii = Ii(z)/f(z) be the total power 
at z divided by the focal spot area, with focusing factor 
/ = A(zf)/A(z) < 1. We typically employ for / the 
result for the on-axis intensity of a gaussian beam [22 : 



/ = [! + (*- zpfjzl 



(53) 



where zq is an effective Raylcigh range. For a Gaussian 
beam with optics F-number F, zq — (A/tt)XF 2 . This 
form approximately fits the random phase plate (RPP) 
smoothed beams designed for NIF (for an appropriate 
zq)- 

Substituting (I , h) = f ■ (Jo, h) into Eqs. ( 47|48 l, and 
freely commuting / with d z , yields the principal equa- 
tions solved by deplete: 



UJ 



d z h{z,ui) = Kiii - Si - 7 (n -hTiH) 



k I ~I I dui — (ti +riii),(54) 



(55) 



r x = ft x and Si = f- x t x . In Eqs. |54p5| and hence- 
forth, all Ii and i\ are understood to have suppressed 
over-tildes, that is, to refer to total transverse powers 
over focal-spot area. Similarly, the plasma-wave ampli- 



tude from Eq. ( 15 1 can be written 

2 



ri2 



2 e 



ck 2 



f a a\" 



(56) 



with of = I i X;/(P em r] i ); see Eq. (|5j. 

All symbols in Eqs. ( 54|55 1 are positive, except Ti may 
be negative for SBS in case ui' 2 < 0. This corresponds to 
the scattered wave having a higher frequency than the 
pump, in the plasma frame. The scattered wave then 
gives energy to the pump, and deplete handles this sit- 
uation correctly. 



III. RAY METHODOLOGY AND MODEL 
LIMITS 

deplete calculates LPI along given plasma conditions 
for a ID profile. A typical application is to study a 
laser beam propagating through conditions given by a 
rad-hydro simulation. We use many independent rays to 
model the whole beam, which introduces some statistical 
inaccuracy. The rays are generally found by tracing 3D 
refracted paths through the rad-hydro output. Although 
strictly not a part of deplete, this is the major way we 
utilize geometric-optics rays. Wave-optics effects, such 
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as laser speckles and diffraction (of both the pump and 
scattered light), are also not included in deplete. We 
present one way to approximate gain enhanacement due 
to speckles in Sec. |Vff| 

However, laser intensity is not found from a rad-hydro 
simulation. Such codes generally treat a laser beam as a 
set of rays, which are absorbed as they trace out refracted 
paths. The laser intensity in a zone is found by divid- 
ing the total power of all rays crossing that zone by its 
transverse area. This approach suffers from several prob- 
lems for our purposes, including the fact that intensities 
remain finite at caustics only due to the finite number of 
rays and zone size. Instead, we run deplete separately 
for each ray, and use a model for the laser beam to give 
an initial intensity (at a sufficiently low density that little 
absorption has occurred) and z-dependent focusing factor 
(generally based on vacuum propagation). The intensity 
along a deplete ID profile is thus independent of refrac- 
tion that occurs due to the plasma. Refractive changes 
in beam intensity occur, for instance, when a beam prop- 
agates between two high-density regions. However, our 
independent-ray treatment has the benefit that caustics 
pose no problem. 

deplete assumes that the laser and scattered light 
follow the same path, and thus see the same plasma con- 
ditions. The two light waves refract differently if their 
wavelengths differ, as in SRS, or in SBS for certain trans- 
verse plasma flows |23j . The departure of ray paths be- 
comes significant when the two rays see sufficiently dif- 
ferent plasma conditions in the gain region for a given 
wavelength that the coupling or other coefficients differ 
significantly. This requires sufficiently strong transverse 
plasma gradients. 



IV. NUMERICAL METHOD 



Let the operator B solve the "bremsstrahlung" system 



We solve the deplete system Eqs. ( 54p5 1 from the 
laser entrance z = to the right edge z — L z . For 
backscatter (considered in this paper), we give Iql and 
*ik(wi) as boundary conditions, where Jl = f(z = 0) 
and fn = f{z = L z ). We solve this two-point bound- 
ary value problem via a shooting method, marching from 
right to left. We guess Iqr and solve the initial value 
problem from z = L z down to z = 0, and iterate un- 
til the resulting Iql is sufficiently close to the desired 
value. Because Iqr is just one scalar, it is more feasible 
to shoot on it than on the set of values iii(wi). General- 
izing our approach to 3D , where one would have to shoot 
on Ior(x,u) over a transverse plane, is much more diffi- 
cult; a different technique for 3D pump depletion is used 
in the code SLIP [23]. For the right-boundary seed value 



d z Io — —kqIoj 
d z i 1 = - Sx, 



(57) 
(58) 



and the operator C solve the "coupling- Thomson" system 



d z Ia — —Iq 
d z ii = -I (ti +T 1 i 1 ) 



dw x — (n +r 1 « 1 ), 

UJi 



(59) 
(60) 



To advance the solution from the discrete gridpoint z n 
down to z™ _1 (the decreasing index matches deplete's 
right-to- left marching), we first apply B for a half-step, 
then C for a full step, then B for a half-step again. The 
splitting theorem guarantees that if B and C are second- 
order accurate operators, then the overall step is second- 
order accurate. Schematically, a complete step is 



{/o.ii}"- 1 =B 1/2 C 1 B 1/2 {I ,i 1 } n . 



(61) 



In usual applications we are given plasam conditions, 
and thus the coefficients in the deplete equations, only 
at a discrete set of points {z n }. We use linear interpo- 
lation to find the coefficients at the needed intermediate 
points, as shown below. We stress that the numerical ac- 
curacy of deplete is strongly influenced by the quality 
of the given plasma conditions. 



A. The bremsstrahlung step B 



B must solve Eqs. ( 57p)8 1 with Ki and Si constant, 
to at least second-order accuracy. This linear system is 
readily solved analytically. Since there are two "half- 
steps" of B in Eq. (61 1, we consider a generic step of size 
Az with initial conditions {Iojii} 1 , yielding new values 
{Io,h}° . X 1 / 2 = (X°+X 1 )/2 denotes the zone-centered 

1/2 

value of some quantity X. If k x ^ 0, we find 



7° = 4exp[ K l /2 Az}, 



i\ = (ij-i 



OT,l/2 



)exp[ 



-k] /2 Az] 



.OT.l/2 



(62) 
(63) 



Eq. ( 63 1 applies separately at each u>\ . For the special 



T72 



0, Eq. (p3| is replaced with 



s} /2 Az 



{n 1 ' 2 = 0). 



(64) 



The rightmost B in Eq. (61 1 advances the system from 
"* to z™ -1 / 2 . Accordingly, for this step, the needed co- 



llR, 



(39 1. The choice seems to have little effect, since volume from z n to z™ _1 : X 1/jl = \(1/4)X 



we either use or the optically-thick from Eq. efficients in Eqs. fl62p4"| are interpolated at 1/4 the way 
' ' ' " ' ' " " " 1 +■ (3/4)X n ] . Simi- 

larly, the leftmost B in Eq. (|6l| ) advances the system from 
z"- 1 / 2 to z™- 1 and uses X^= [(3/4)!"- 1 + (l/4)X n ]. 
In both cases Az = (z™ - z"- 1 )/2. 



sources (either TS or bremsstrahlung) typically produce 
a comparable or larger noise level after a short distance. 
We solve Eqs. ( 54p5 l by operator splitting [251 125] - 
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B. The coupling-Thomson step C 

We now turn to the C operator. Iq is evolved via a 
conservation law of the C system, Eqs. ( 59||60 1: 



, w . 



(65) 



where i\l is the solution to the full deplete equations. 
Gi and Gd are exactly equal under the following condi- 
tions: there is no pump depiction, no TS (ti = 0), no 
absorption of scattered light {n\ — 0), and no volume 
bremsstrahlung noise (Si = 0); the only seeding in de- 
plete is then via the boundary values iir{u>i). 



On the discrete z grid, this gives 



SRS benchmark 



rn- 1 jn 



dwi — (*r 



(66) 



Before doing this, we must advance i\ using Eq. (60) 



with constant Iq = Iff (that is, we neglect pump deple- 
tion within a zone). This gives rise to a numerical chal- 
lenge. Namely, the coefficients T\ and Ti are both pro- 
portional to |e| -2 , and contain a narrow resonance where 
Re e = if Im e is small (that is, where the beating of the 
light waves drives a natural plasma wave). Integrating 
through these sharp peaks with a standard ODE method 
like Runge-Kutta performs very poorly unless the res- 
onance is well-resolved by the z grid (which it usually 
is not). To alleviate this problem, the key observation 
is that e itself varies slowly in space, even though |e|~ 2 
varies rapidly near resonance. We can therefore represent 
e as linearly varying with z across a cell, and analytically 
solve the resulting system. We merely quote the result 
here, and refer the reader to Appendix B for the deriva- 
tion and definition of the relevant quantities: 



(*i +V)e 



BrAw n 



(67) 



V. BENCHMARK ON LINEAR PROFILES 

This section compares the results of deplete with 
those of NEWLIP and pF3D on two contrived profiles with 
weak linear gradients, one for SRS and another for SBS. 
deplete and pF3D embody quite different physical mod- 
els, each with their own approximations and limitations. 
One can view their favorable comparison here as a "cross- 
validation" of these models in a regime where they should 
agree. 

To compare with the NEWLIP linear gain G\ (see Ap- 
pendix A), we need a noise level against which to compare 
the deplete scattered spectrum at the laser entrance, 
i\L- For this noise level we choose i\ r at z = 0, given 



0): 



by solving Eq. (551 with just the bremsstrahlung terms 



■br 
Kill 



Si. 



(68) 



This is exactly Eq. ( 58 ) . We then introduce the deplete 
gain G d : 



1 1L 



"scattering'' 
"noise" 



(69) 



The spatial profiles of our SRS benchmark plasma con- 
ditions are shown in Fig. [T] We use a profile length 
L z = 510Ao, pump vacuum wavelength Ao = (1054/3) 
nm, fully-ionized H ions with Ti — I keV, and no plasma 
flow (u = 0) . In both the deplete and pF3D runs of this 
section, SBS was not included. Fig. [2] plots the result- 
ing reflectivities for several pump strengths. Although 
these are all above the homogeneous absolute instabil- 
ity threshold of Ig b « 0.21 PW/cm 2 , the time-dependent 
pF3D runs rapidly approach a steady state and show no 
signs of a temporally-growing mode [42] . The weak gra- 
dients, or incoherent noise source, may lead to stabiliza- 
tion. After increasing exponentially with Iql for weak 
pumps, the reflectivity rolls over. This saturation due to 
pump depletion is generic for three-wave interactions in 
the strong damping limit, as demonstrated analytically 
by Tang [27]. 

We compare the gains Gi and Gd from NEWLIP and de- 
plete, for several pump strengths, in Fig. [3} The general 
shapes of the gains are quite close, although their abso- 
lute levels differ. For the weakest pump strength, where 
pump depletion plays little role (as can be inferred from 
the reflectivity plot in Fig. [2j, the peak Gd is slightly 
higher than G;. This is due to the volume sources in de- 
plete, namely TS and bremsstrahlung noise. To illus- 
trate this, we plot Gd found with no Thomson scattering 
(ti = 0) as the black dotted curve. It lies between the 
two other curves near the peak, and overlaps Gi away 
from the peak. The curves for the two larger values of 
Iql in Fig. [3] show Gd to be progressively farther below 
Gi at peak. This results from pump depletion, which the 
reflectivity plot clearly shows is significant for Iql > 0.8 
PW/cm 2 . The bremsstrahlung noise level i\ r varies be- 
tween (2.4-4.1) xl0~ 9 W/cm 2 /(rad/sec) over Ai = 650 
to 550 nm. 

We also compared deplete to the massively-parallel, 
paraxial laser propagation code pF3D [25]. This code 
solves for the slowly- varying envelopes of the pump laser, 
nearly-backscattered SRS and SBS light waves, and the 
daughter plasma waves, in space and time. A carrier 
Lo en is chosen for each mode (except for the ion acoustic 
wave), and the corresponding rapid time variations are 
averaged over. A local eikonal k en , given by the appro- 
priate uj en and dispersion relation with local plasma con- 
ditions, contains the rapid space variation. Kinetic quan- 
tities, such as Landau damping rates and Thomson cross- 
sections, are variously found from (linear) kinetic formu- 
las or fluid approximations. There is no bremsstrahlung 
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FIG. 1: Plasma conditions for SRS benchmark. 



hq+o _i i i i I i i i i i i i i i I i i i i i i i i i L. 




— i ■ i ■ I • i ■ i ■ i • i • I ■ i ■ i • i ■ i ■ | — 

0-5 1.0 1.5 



Iol (PW/cm 2 ) 

FIG. 2: (Color online.) SRS reflectivity vs. pump intensity 
for the SRS benchmark profile of Fig. [l] The black circles 
and red squares are for pF3D and deplete, respectively. 



source, but the pump and scattered light waves all ex- 
perience inverse-bremsstrahlung damping. The plasma 
waves undergo Landau damping, and the advection term 
v g 2d x ri2 is retained (i.e., they are not treated in the 
strong damping limit). The noise source in pF3D is 
plasma-wave fluctuations chosen to produce the correct 
TS level, and uniformly distributed over a square in k± 
space (corresponding to the transverse x and y direc- 
tions) extending to half the Nyquist k in both k x and 
k y . 

To replicate the ID model of deplete, we performed 
"plane- wave" simulations in pF3D. The incident laser 
at the z = entrance plane is uniform in the x and y 
directions (i.e., there is no structure like speckles), both 
of which are periodic with size L x = L y = 128Ao and grid 
spacing dx = dy = 1.33Ao- The z spacing is dz — 2\q. 
As described above, the TS noise fills a square in k± 
space extending to k x ,k y — ±fci n , with k\ n — (3/16)fcot> 
and k 0v = uo /c. We enveloped the SRS backscattered 
light around ujf n = 0.592wo (Ai=593.3 nm), which has 
the highest linear gain. Over the slight variation of our 
profile, the average kf n = 0A61k Ov . 

deplete requires a solid angle il c , which we express 
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FIG. 3: (Color online.) deplete gain Gd (black solid), 
NEWLIP gain Gi (red dashed), and Gd with no TS for Iol = 0.4 
PW/cm 2 (n = 0, black dots), for SRS benchmark. TS 
and volume bremsstrahlung noise enhance Gd over Gi for 
the smallest Iol, while pump depletion suppresses Gd for the 
larger two. 

in terms of an F-number F, for TS and bremsstrahlung 
emission (we excluded the latter for pF3D comparisons). 
Taking kf n and k\ n to determine the focal length and 
spot radius, one finds F — kf n /2k\ n — 1.23. The scat- 
tered light does not uniformly fill the noise square in k± 
space, but rather develops into a somewhat hollow "ring" 
with a radius ~ 0.12fcoi, (departing more from a square 
for stronger pumps); there is some ambiguity in the ap- 
propriate F to use. We choose F = 1, which leads to 
very close reflectivities for the weakest-pump case shown 
in Fig. [2] and is near the noise-square estimate F = 1.23. 
Sidescatter at these angles may stress the accuracy of 
Pf3d's paraxial approximation. 

Figure [2] shows the deplete and pF3D SRS reflectiv- 
ities for the benchmark profile. The pF3D values are 
taken at t =39.4 ps, after which time all reflectivities 
remain roughly constant (the laser ramped from zero to 
full strength over 10 ps). The agreement is quite good, 
especially in the linear (weak pump) and the strongly- 
depleted (strong pump) regimes. This increases confi- 
dence in the validity of the different approximations made 
in both codes. It took about 2 sees of wall time for de- 
plete to run on one Itanium CPU, as opposed to 5300 
sees on 16 of these CPUs for pF3D to advance 10 ps. 



B. SBS benchmark 

We performed an SBS benchmark (with SRS ne- 
glected) using the profiles in Fig. [2] The ions were fully- 
ionized He (Z = 2, A = 4) with T, = T e /5. The parallel 
flow velocity u is shown normalized to the local acous- 
tic speed (? a = {ZT e + 3Ti)/Am p . The pump wavelength 
and profile length match the SRS benchmark. The SBS 
reflectivity vs. pump strength is plotted in Fig. [5] which 
shows pump depletion for Iol ?J 1-25 PW/cm 2 . We es- 
timate the absolute threshold lQ b — 2.6 PW/cm 2 and 
stay below this. We used F = 1.7 since this gives good 
agreement with pF3D "plane-wave" simulations for low 



9 




z / Xo z / Xo 



FIG. 6: (Color online.) SBS deplete gain Gd (black solid), 
FIG 4- SBS benchmark profile newlip gain G t (blue dashed) , and G d without TS for 

Iol =0.6 PW/cm 2 (n = 0, black dotted), for SBS bench- 
mark profile. 
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FIG. 5: SBS reflectivity for SBS benchmark profile. The 
squares are deplete results, and the dashed line is an ex- 
tension of the low- Jql results. 



However, for larger values of Jo a ring in k± space 
develops, similar to the SRS runs, and is accompanied by 
a large increase in reflectivity. 

Figure [6] compares the deplete and newlip gains, Gd 
and G/. For the smaller two pumps we see the enhance- 
ment of Gd over Gi due to TS (even though pump deple- 
tion has set in for the second case Iol — 1-4 PW/cm 2 ), 
as discussed in Sec.|Vl] The dotted black curve for Iol = 
0.6 PW/cm 2 is Gd computed with no TS, and shows the 
modest increase in Gd stemming from bremsstrahlung 
volume (as opposed to boundary) noise. The elevated 
plateau of Gd to the left of the peak is also due to 
TS. Iol =2.5 PW/cm 2 gives G d < G t due to strong 
pump depletion. In all cases the wavelength and width 
of the main peak of the two spectra are similar. i\ r , 
the bremsstrahlung solution, varies slightly from (4.17- 
4.25) xlCT 9 W/cm 2 /(rad/sec) over A x - A = 20 to -3 

A. 



VI. THE RELATION OF THOMSON 
SCATTERING TO LINEAR GAIN 

As seen in our benchmark runs, TS leads to an en- 
hancement of the deplete gain compared to the newlip 
gain (for negligible pump depletion) . This is readily seen 
via the scattered-wave equation with just coupling and 
TS, Eq. (pbl: 



d z h = -loin + Ftix). 



We use Eq. (51 1 to obtain 



d z ii = ~7(*t + h)- 



(70) 



(71) 



7 = Io^i is the spatial gain rate. Typically, 7 has a 
narrow peak in z at the resonance point, while i T varies 
slowly. For simplicity, we hold i T constant at the reso- 
nance point, and solve for i\ across the region z = to 
L, which includes the resonance. In our usual notation, 



l\L 



(iiR + v)e Gi - i T . 



(72) 



Gi = f Q z dzj is the newlip linear gain. For Gi 1, 
i\L = iiu(l + Gi) + i T Gi, and emission due to the bound- 
ary source dominates over TS. In the opposite limit, 



ilL = (iiR + ir)e Gl , 



e Gl > 1. 



(73) 



TS therefore gives rise to an effective boundary source 
i T (for a narrow resonance). In this sense, it does not sig- 
nificantly alter the shape of the gain spectrum (i r varies 
slowly with loi). However, it does lead to a difference 
in the absolute magnitude of the scattered spectrum, as 
embodied in an "absolutely-calibrated" gain like Gd- As 
an illustration, let us take i\R — i® T , the optically-thick 
bremsstrahlung result of Eq. ( 39 ) , for simplicity evalu- 
ated at the resonance point in the Jeans limit huji <^ T e . 
Moreover, we set Ti = T e so that i T assumes the simple 



form of Eq. (52 1. The effective seed is then 



ii R + i r 



■OT 



^1 , p Wo U>o 



V. 



UJi LU 2 



(74) 
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The second term on the right (= i T /ip T ) is typically < 10 
for SRS: for our SRS benchmark, i r /if T w 3. But, it can 
be quite large for SBS since ojq ^S> u>' 2 (for our SBS bench- 
mark, i T /ii « 400). A similar result is found in Ref. [BJ. 
The authors explain this on the thermodynamic ground 
that bremsstrahlung and Cerenkov emission (which pro- 
duces TS) generate equal light- and plasma-wave action, 
so the light-wave energy dominates by the frequency ra- 
tio. This manifests itself in the u)o/u>' 2 factor in Eq. (74 1, 
which is much larger for SBS. 



VII. SIMULATION OF SBS EXPERIMENTS 

Experiments have been conducted recently at the 
OMEGA laser to study LPI in conditions similar to 
those anticipated at NIF [25]. These shots use a gas- 
filled hohlraum, and a set of "heater" beams to pre- 
form the plasma environment. An "interaction" beam is 
propagated down the hohlraum axis after being focused 
through a continuous phase plate (CPP) [3D] with an 
f/6.7 lens to a vacuum best focus of 150 /im. The plasma 
conditions along the interaction beam path have been 
measured using TS [3T], validating 2-dimensional hydra 
[32J hydrodynamic simulations that show, 700 ps after 
the rise of the heater beams, a uniform 1.5-mm plasma 
with an electron temperature of s»2.7 keV [55] , 

Figure f7fa) displays the instantaneous SBS reflectivity 
increasing exponentially with the interaction beam in- 
tensity 700 ps after the rise of the heater beams. These 
experiments employed a 1 atmosphere gas-fill with 30% 
CH 4 and 70% C3H8 to produce an electron density along 
the interaction beam path of 0.06n c o. Three-dimensional 
pF3D simulations agree well with the experiments [51] . 
Unlike the plane- wave simulations discussed in Sec. |V A| 
these simulations include the full speckle physics. The 
deplete results (blue solid curve) fall well below the 
experimental data in the regime where pump depletion 
does not play a significant role (Iq < 2 PW/cm 2 ). This 
indicates that speckles are enhancing the SBS. 

One way to approximate the speckle enhancement is 
to consider how much the coupling increases for the 
completely phase-conjugated mode [35]. This mode has 
a transverse intensity pattern perfectly correlated with 
that of the pump, over several axial ranks of speckles, 
and therefore enhances the coupling coefficient Ti [55] . 
For an RPP-smoothed beam with intensity distribution 
~ e _I / Jt! , this effectively doubles r x . This should pro- 
vide an upper bound on the reflectivity so long as the 
gain per speckle is < 1. If this is not the case, the gain in 
a speckled pump suffers a mathematical divergence (mit- 
igated by pump depletion) as described in Ref. |37j . Our 
phase-conjugate considerations would then not apply. 

The blue dashed curve in Fig.[7]shows the deplete re- 
sults with twice the nominal coupling. The 2 x Ti curve is 
always above the experimental reflectivities. The thresh- 
old intensity for which SBS equals 5% is 1.8 PW/cm 2 and 
0.9 PW/cm 2 for deplete with the nominal and twice- 
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FIG. 7: (Color online.) (a) SBS reflectivity for OMEGA ex- 
periments with CH gas fill and T e w 2.7 keV (described in 
text). Black circles are measured values, the blue solid curve 
is deplete calculations with the nominal coupling Fi, and 
the blue dashed curve is deplete calculations with 2 x fi. 
(b) deplete and pf3d SBS reflectivities for a similar config- 
uration but T e w 3.3 keV. Black crosses are pf3d simulations, 
and the blue curves are the deplete results as in (a). 



nominal coupling, respectively, while the experimental 
threshold is wl.5 PW/cm 2 . 

Comparison of deplete and pF3D is displayed in Fig. 
|7][b). These calculations were performed using plasma 
conditions from a hydra simulation, for a configuration 
similar to that of Fig.[7]ja), but with a higher heater-beam 
energy. The resulting conditions are similar, except the 
electron temperature is higher (about 3.3 keV). The de- 
plete reflectivity with the nominal coupling (solid blue 
curve) lies below the pF3D results for the two interme- 
diate values of Iq. This demonstrates speckle effects en- 
hance the pF3D reflectivity for moderate Iq. The de- 
plete results for 2 x Ti (dashed blue curve) are always 
above the pF3D results. Preliminary analyses with de- 
plete and pF3D of OMEGA experiments designed to 
study ion Landau damping in SBS [ TO] also show a sig- 
nificant enhancement due to speckles. 
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FIG. 8: Nominal intensity at best focus for 285 eV NIF ig- 
nition design ("NIF example"), found by dividing the laser 
power per quad by the focal spot size. The peak intensity 
corresponds to 6.9 TW/quad. 




I i l i i ■ l ■ l i I i l i l 



10 



285 eV Point Design 



230 



refraction 



| M i I i I i J < | i I ■ I i I i I ■ | i m- 

-0.5 00 0.5 

z (cm) 

FIG. 9: (Color online.) Materials and laser beam cones for 
NIF example. 



VIII. ANALYSIS OF NIF IGNITION DESIGN 

In this section, we exercise deplete on an actual NIF 
indirect-drive ignition target design. The target was de- 
signed using the hydrodynamic code lasnex [38 . For 
more details about the design see Ref. [39]; LPI analysis 
for this and similar ignition targets, including massively- 
parallel, 3D pF3D simulations, can be found in Ref. [30] , 
The design utilizes all 192 NIF beams (at 351 nm "blue" 
light), which deliver 1.3 MJ of laser energy. We analyze 
LPI along the 30° cone of beams (one of the two "in- 
ner" cones). The pulse shape for one quad (a bundle of 
four beams), expressed as nominal intensity at best fo- 
cus, is shown in Fig. [8] and reaches a maximum of 0.33 
PW/cm 2 . The speckle pattern for a quad approximately 
corresponds to an F- number of F = 8, which we use for 
deplete's noise sources (but each beam individually has 
F = 20 optics) . The focal spot is elliptical with semi-axis 
lengths of 693, 968 fj,m. The peak temperature of the ra- 
diation drive is 285 eV. The materials are as follows: the 
capsule ablator is Be, a plastic (CH) liner surrounds the 
laser entrance hole, the hohlraum wall is Au-U with a 
thin outer layer of 80% Au-20% B (atomic ratio), and 
the initial fill gas is 80% H-20% He. The lower-Z com- 
ponents are included in the last two mixtures to reduce 
SBS by increasing the ion Landau damping of the acous- 
tic wave. 

We performed deplete calculations, with both SRS 
and SBS, at several times and over 381 ray paths for 
each time. One must take an appropriate "average" over 
the rays to characterize the LPI on a cone. Regarding 
newlip gains, this has led to several approaches. These 
include averaging the gain, finding the maximum gain, 
or averaging e Gl . This last method stems from assuming 
there is no pump depletion and noise sources are inde- 
pendent of scattered frequency; in this limit, the reflec- 
tivity should be roughly proportional to e G ' . However, 
this averaging, and a fortiori taking the maximum, can 



be dominated by gains that are larger than physically al- 
lowed by pump depletion or other nonlinearity. One can 
attempt to include pump depletion via a Tang formula 
for Gi at each w\ [37] . 

deplete allows for more physical ray-averaging 
schemes. To the extent the transverse intensity pattern 
of a cone is uniform, each ray represents the same inci- 
dent laser power. Averaging deplete's ray reflectivities 
then measures the fraction of incident power that gets 
reflected. Pump depletion is of course included, which 
limits backscatter along high-gain rays in a physical way. 
The reflectivities and scattered spectra plotted here are 
simple averages over the rays. 

The reflectivities for several times near peak laser 
power, for the 30° cone, are shown in Fig 
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The re- 
sults for three different cases are presented. First, the 
solid lines give the reflectivities computed with the un- 
modified deplete equations. To quantify the role of 
re-absorption of scattered light in the target, we re-ran 
deplete with m = 0. This leads to the dashed lines. 
Finally, to bound the enhancement due to speckles, we 
plot the results when T\ is doubled (and K\ ^ 0) as the 
dotted lines. 

The spectra of escaping SRS and SBS light (averaged 
over rays) are shown in Fig. pTp^2| The SBS feature at 
a wavelength shift of 5-8 A comes from the Be ablator 
blowoff. A much weaker feature appears from 12-13 ns 
at 12-15 A, and occurs in the gas fill. The SRS spectrum 
is more irregular, showing two main features separated 
by w20 nm that move to higher Ai as time increases. 
In addition, there are narrow features at higher Ai that 
originate near the hohlraum wall; these would be reduced 
in a ray-averaged gain, since the exact Ai active for each 
ray depends sensitively on conditions near the wall and 
therefore varies from ray to ray. Re-absorption strongly 
suppresses these high- A i spikes, as is seen in the SRS 
spectra with and without re-absorption at t — 13.75 ns 
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FIG. 10: (Color online.) deplete SRS and SBS ray-averaged 
reflectivities I\l for NIF example. Solid lines are the nominal 
case (re-absorption and Fi unsealed), dashed lines are the 
nominal Ti but no re-absorption of scattered light («i = 0), 
and dotted lines are 2 x Ti with re-absorption. 
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FIG. 12: (Color online.) SBS streaked spectrum for NIF ex- 
ample, nominal case (ki 7^ 0, 1 x Ti). The white-yellow streak 
from 5-8 A occurs in the Be ablator, while the weaker feature 
from 12-15 A occurs in the gas fill. 
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FIG. 13: (Color online.) deplete SRS spectrum at time 
13.75 ns for NIF example, smoothed over ~ 1 nm. The black 
solid and red dashed lines are computed with («i 7^ 0) and 
without (ki = 0) re-absorption of scattered light, respectively. 



FIG. 11: (Color online.) SRS streaked spectrum i\t for NIF 
example, nominal case (k\ / 0, 1 x Ti). 



in Fig. 13 Collisional plasma-wave damping, currently 
not in deplete, may reduce the high-Ai scattering (the 
Landau damping of the \0w-k2\De plasma waves is neg- 
ligible). 

Besides backscatter, deplete also provides the pump 
intensity Io(z) along each ray. This indicates how much 
laser energy is transmitted to a given location, which is 
a crucial aspect of a whether LPI degrades target per- 
formance. In cases where the backscattered light under- 



goes significant absorption as it propagates out of the 
target (as happens to SRS for the design analyzed here), 
the measured reflectivity can understate the level of LPI. 
The laser transmission can reveal this fact. Figure 14 'a) 



presents Io, averaged over all the rays, at a given n e . 
This is a ID presentation of how much energy reaches a 
given density, although in the full 3D geometry different 
rays reach the same n e at different locations. I with 
just pump absorption, as well as the deplete solutions 
with pump depletion for the nominal case and 2 x Li, are 
shown. Pump depletion is barely discernible in the nomi- 
nal case, but is significant in the 2 x Ti case. For instance, 
in the latter case Iq at n e /n C Q = 0.2 is only 60% of its 
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FIG. 14: (Color online.) (a) Laser transmission for NIF exam- 
ple at 12.5 ns (peak power): black solid curve is the nominal 
deplete solution with pump depletion, red dashed curve is 
with just inverse-bremsstrahlung absorption, and black dot- 
ted curve is the deplete solution with 2 x Fi. (b) SBS (blue) 
and SRS (red) scattered intensities for the nominal deplete 
solution. Calculation of intensity at a given n e is described 
in text. 
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FIG. 15: (Color online.) SRS spectral density i\ vs. n e /n c o 
and Ai, in decibels, at 12.5 ns (peak power), for NIF example. 



absorption-only value. The wavelength-integrated SRS 
and SBS I\ are shown in Fig. 14 b), and the scattered 
spectra vs. n e are shown in Figs. |15|l6l SRS in particu- 
lar develops at several different densities, corresponding 
to different wavelengths, as can be seen in Figs. [13] and 



20 J 



< 10- 



0- 



1 

00 



0.1 



0.2 



-100 



ED 



FIG. 16: (Color online.) SBS spectral density i\ vs. n e /n c o 
and Ai — Ao, in decibels, at 12.5 ns (peak power), for NIF 
example. 



IX. CONCLUSIONS AND FUTURE 
PROSPECTS 

We have derived a ID, steady-state, kinetic model for 
Brillouin and Raman backscatter, that includes pump de- 
pletion, bremsstrahlung damping and fluctuations, and 
Thomson scattering. This model is implemented by the 
code deplete, which we have presented as well. This 
work extends linear gain calculations, by including more 
physics while retaining its low computational cost. In 
particular, deplete provides the scattered-light spec- 
trum and intensity developing from physical noise, which 
can be compared against more sophisticated codes and 
experiments. The transmitted pump laser along the pro- 
file is also found, which is important for assessing an ICF 
target design, especially when re-absorption of scattered 
light reduces the escaping backscatter from its internal 
level. 

We presented benchmarks of deplete on contrived, 
linear profiles, as well as analysis of OMEGA experi- 
ments and a NIF ignition design. The benchmarks re- 
veal the deficiencies of linear gain, namely the neglect 
of TS, pump depletion, and re-absorption. Comparisons 
with pF3D provide a cross-validation of the two codes in 
a regime where they should agree. The OMEGA SBS 
experimental data, as well as pF3D simulations of these 
shots, show much more reflectivity than deplete gives, 
for intensities where pump depletion is weak. This en- 
hancement is due to speckle effects. We showed an upper 
bound on this enhancement is given by doubling the de- 
plete coupling coefficient T\, which comes from consid- 
ering the phase-conjugated mode in an RPP-smoothed 
beam. The ignition design analysis gives reasonably low 
backscatter levels for the nominal laser intensity and in- 
cluding re-absorption, with SRS dominating SBS. How- 
ever, if re-absorption is neglected, or especially if Ti is 
doubled, the backscatter appears more worrisome. The 
laser transmission supports these conclusions. 
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Ray-based gain calculations have been used for some 
time to model LPI experiments, and deplete can pro- 
vide more detailed comparisons. An early application of 
gain to hohlraum targets is Ref. [H], where hohlraums 
filled with CH gas were driven by laser beams with and 
without PS and SSD. Without SSD, reasonable agree- 
ment was found between measurments and the time- 
dependent SBS gain spectrum. However, there was a 
large difference in peak SRS wavelength between mea- 
surements and the gain spectrum, which may be due to 
laser filamentation changing the location of peak SRS 
growth. 

Several future directions exist for deplete. One is 
to include an "independent speckle" model for gain en- 
hancement, where one solves the deplete equations over 
a speckle length for a distribution of pump intensities and 
then re-distributes the power. This would not describe 
correlations among axial ranks of speckles, caused e.g. by 
phase conjugation, deplete also enables some new diag- 
nostics and applications. The pump and scattered inten- 
sities found by deplete can be used to compute the local 
material heating rate due to absorption. This could be 
incorporated into a hydrodynamic code, thereby coupling 
LPI to target evolution in a self-consistent, if simplified, 
way. In addition, the plasma- wave amplitudes found by 
deplete can be compared against thresholds for various 
nonlinearities to assess their relevance, and may allow 
estimation of hot electron production by SRS. 

Despite its promise, there are limits inherent to any 
ID or ray-based approach, stemming from 3D wave op- 
tics (e.g. diffraction, speckles, filamentation, and beam 
bending) . A 3D paraxial code called slip [24 , which like 
deplete operates in steady state and uses kinetic coef- 
ficients, is being developed. This model is in some sense 
intermediate between deplete and pF3D. ID codes 
like deplete still have a valuable role. They can ana- 
lyze hundreds of rays, using hundreds of scattered wave- 
lengths, in ~ minutes, thus allowing designs to be rapidly 
analyzed and compared. The resulting time-dependent 
spectra allow for contact with experimental diagnostics, 
and are frequently needed, for example, to choose the 
carrier k and to for pF3D. 

Laser-plasma interactions have proven to be a very 
challenging area of plasma physics, owing to the variety 
of relevant physics and extreme range of scales involved. 
This has led to an equally extreme range of modeling 
tools, from ID gain estimates to 3D kinetic simulations. 
By fully exploiting these tools, each with their uses and 
limitations, a more complete picture is emerging. 
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APPENDIX A: NEWLIP 

In this appendix, we document the laser interaction 
post-processor newlip, of which deplete can be viewed 
as an extension. The equations underlying newlip are 



d z I (z) 
d z i\(z,ui) 



(Al) 
(A2) 



The first of these is Eq. (54 1 with no pump depletion 
(ti = = 0), and the second is Eq. (551 with no 
bremsstrahlung or TS (ki = Ex = t\ = 0y That is, 
only the absorption of the pump, and coherent coupling 
to scattered light waves, are modeled. The boundary con- 
ditions are Iq(z = 0) = Iql (the known pump at the laser 
entrance), and i\(z — L Zi uj\) — 1. We thus solve for a 
unit scattered- wave boundary seed, which is permissible 
for this linear system. 

We readily solve Eqs. ( [Al|A2| to find 



Hz) 
h(z) 

G t (z) 



- p Gi(z) 



dz 1 Tx{z>)I Q (z>). 



(A3) 
(A4) 

(A5) 



Gi(z) is the linear intensity gain exponent, and is 
newlip 's main result. The total gain across the profile 
isG ; (z = 0). 

The numerical computation of Gi suffers from the 
problem of narrow resonances, similar to deplete. The 
coupling coefficient Ti (see Eq. ( fl~7| ) is sharply peaked 
near the resonance point where Ree = 0. newlip ad- 
dresses this challenge in a way analogous to how deplete 
handles the coupling-Thomson step, as outlined in Ap- 
pendix B. In particular, the integration of Eq. (A2) from 
z n down to can be cast in the form 



In- 



= ImSo, 



s 



(A6) 



with S = -Io/r s Xe(l + Xi)- Although S(z)/e(z) is 
sharply-peaked near resonance, S(z) and e(z) are them- 
selves slowly-varying with z. We approximate S(z) 
S n-i/2 + ( z _ z"- 1 /2)AS*/Az (and similarly for e), with 
X n-i/2 = ( Xn + x"-!)/2 and Al = (X n - A"" 1 ) 
for some quantity X . With this representation, and 
X = X n ~V 2 /AX, we find 



Sq 



AzAS 
Ae 



1 + (e - S) In - 



(A7) 



This formula is valid provided either |Ree| > 1/2 or 
Ime 7^ 0. For accuracy, we also want |Ae| to not be too 
small (which obtains, e.g., for a flat profile). We there- 
fore require |e| to be less than some large number. If 
any of these conditions does not hold, we simply assume 
S = S n - X / 2 and e = e"" 1 / 2 across the cell to find 



Sn = A, 



' ,n-l/2 ' 



(A8) 
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APPENDIX B: NUMERICAL SOLUTION OF 
THE COUPLING-THOMSON STEP 



This appendix provides a derivation of Eq. (67), the 



The linear change of variable 



solution f or i \ in the coupling- Thomson step. We must 
solve Eq. (60 1, from z n down to z™" 1 , with Iq = Ift and 
all coefficients except |e| 2 evaluated at z n_1 ' 2 . We write 
this equation as 



(z - z ) 



(B8) 



transforms Eq. ( Bl I to 



K T + K T i x 

m n-1/2 n-1/2 



K T = I» T 'g -ff- 



AY 



(Bl) 

//«rr 1/2 pr 1/ (fe2) 



d s ii = 
B T r = 



_B T + B T i x 
1 + s 2 ' 

|Im [e^-VZAe*] 



(B9) 
(BIO) 



As mentioned above, the principal numerical difficulty is 
that |e 



is sharply peaked near resonance (Ree = 0). A second chan S e of vari able to w = atan s yields 



Since Re e generally passes through zero slowly, we Taylor 
expand e within each zone and solve the resulting system 
analytically. 

Define the zonal average and difference A' 1 ^ 1 / 2 = 
(1/2)(A™ + A™" 1 ) and AA = A" - A"" 1 for the quan- 
tity A. We expand e about the zone center z™^ 1 / 2 to 
find 



d w ii = -{B T + B r ix). 



(Bll) 



This equation is solved to give the result used in Eq. ( 67 1 

i'l- 1 = (% + i T )e BrAw 



-1/2 



+ zAe 



z = 



z - z"- 1 / 2 
Az 



We can then write 

|e| 2 = e 1 + |A e | 2 (z-z ) 2 , 



ei 



lAeMlm 



£ n-l/2 Ae , 



5 = -|Aer 2 Re e n ~ 1/2 A 



(B3) 
(B4) 

(B5) 
(B6) 
(B7) 



Aw,, = atan 



atan s 



(B12) 
(B13) 



i T = B T /Br is also given by Eq. (51 1 



If Ae is sufficient ly sm all (or zero, as for a flat profile), 
we do not use Eq. p3|, but instead e w e"" 1 / 2 . We can 



then immediately solve Eq. (Bl I to find 



= (»i + V)cxp [A r |e"- 1/2 r 2 Azl - i T . (B14) 
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